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^ I A search for physics beyond the Standard Model in neutral current deep inelastic scat- 

tering at high negative four-momentum transfer squared Q^ is performed in e^p collisions 
at HERA. The differential cross section da/dQ^ , measured using the full HI data sample 
corresponding to an integrated luminosity of 446 pb~^, is compared to the Standard Model 
prediction. No significant deviation is observed. Limits on various models predicting new 
phenomena at high Q^ are derived. For general four-fermion eeqq contact interaction mod- 

5^ I els, lower hmits on the compositeness scale A are set in the range 3.6 TeV to 7.2 TeV. Lep- 

toquarks with masses Mlq and couplings A are constrained to Mlq/A > 0.41 — 1.86 TeV 
and limits on squarks in i?-parity violating supersymmetric models are derived. A lower 
limit on the gravitational scale in (4H-n) dimensions of Ms > 0.9 TeV is established for 
low-scale quantum gravity effects in models with large extra dimensions. For the light 
quark radius an upper bound of Rq < 0.65 • 10^^^ m is determined. 
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1 Introduction 

Deep inelastic neutral current (NC) scattering e^p — >■ e^X at high negative four-momentum 
transfer squared Q^ allows the structure of eg interactions to be probed at short distances and 
to search for new phenomena beyond the Standard Model (SM). Using the concept of four- 
fermion contact interactions (CI) the interference of the photon and Z-boson fields with any 
new particle field associated to larger scales can be investigated. 

Results from searches for contact interactions in ep interactions at HERA have been pre- 
viously reported by the HI [IlllH and ZEUS |l3l|4l collaborations. Therein, genuine contact 
interaction models, models with leptoquarks and supersymmetric scalar quarks (squarks), low- 
scale quantum gravity models with large extra dimensions and compositeness models of quarks 
have been investigated by searching for deviations from the SM expectation at high Q^. Contact 
interaction studies have been also performed at LEP Q. 

Such models have also been investigated in direct searches at HERA, the Tevatron and the 
LHC. Searches for leptoquarks involving lepton flavour violation [6 J and squarks in i?-parity 
violating (f^p) supersymmetric models [7] have been published by the HI collaboration using the 
full HERA data. Searches for leptoquark pair production in proton-proton collisions at a centre- 
of-mass energy of -\/s = 7 TeV were reported by the ATLAS [SJ and CMS [9] collaborations, 
excluding first generation scalar leptoquarks up to 376 TeV and 384 TeV, respectively. These 
results surpass limits obtained at the Tevatron [[TOll . Stringent limits on low-scale quantum 
gravity models with large extra dimensions using di-jet, di-electron and di-photon final states 
have been reported by the D0 collaboration ifTTlfTll and recently, by the CMS collaboration, 
excluding mass scales below 1.6-2.3 TeV depending on the model [fT3l . 

The analysis in this paper is based on the full HI data sample collected in the years 1994- 
2007, which corresponds to an integrated luminosity of £ = 446 pb^^ and represents a factor of 
3(12) increase in statistics for e^p (e~p) collisions with respect to the previous publications flU. 
The same method is used as in previous analysis which is superseded by the results presented 
in this paper. 

2 Contact Interaction Models 

New physics phenomena in fermion-fermion scattering experiments may manifest themselves 
in deviations of the differential cross section do/dQ^ from the SM expectation, and may be 
related to new heavy particles with masses Mx much larger than the electroweak scale. In the 
low energy limit a/s ^ Mx such phenomena can be described by an effective four-fermion CI 
model. Different implementations of this effective model are summarised in the following. 

2.1 General contact interactions and compositeness 

In ep scattering, the most general chiral invariant Lagrangian for neutral current vector-like 
four-fermion contact interactions can be written in the form [|14llT5l : 

^v = ^ Y^ Vlbi^al^,ea){qbl^qb) , (1) 

q a,b = L,R 



where 77^^ are the CI coupling coefficients, a and b indicate the left-handed and right-handed 
fermion helicities and the first sum is over all quark flavours. In the kinematic region of interest 
mainly the valence quarks (u and d) contribute. 

In the case of general models of fermion compositeness or substructure the CI coupling 

coefficients are defined as: 

Air 

Vlb = 4bj^- (2) 

New physics models are then characterised by a common compositeness scale A and the coef- 
ficients e^^, which describe the chiral structure of the coupling and may take the values ±1 or 
0, depending on the scenario, for example pure left-handed (L), right-handed (R), or vector (V) 
and axial-vector (A) couplings. Depending on the model and the sign of the coefficients, the 
new physics processes interfere either constructively or destructively with the SM processes. 



2.2 Leptoquarks 

Leptoquarks, colour triplet scalar or vector bosons carrying lepton and baryon number, appear 

naturally in extensions of the SM which aim to unify the lepton and quark sectors. For lepto- 

quark masses Mlq much larger than the probing scale Mlq ^ ^/s, the coupling A is related to 

the CI coupling coefficients via: 

A^ 
Vlb = 4b^- (3) 

The classification of the leptoquarks follows the Buchmiiller-Riickl-Wyler (BRW) model [fT6l . 
in which the coefficients e^^ depend on the leptoquark type ifTTl and take values 0, ±|, ±1, ±2. 

Two leptoquark types, Sq and 5^2, have quantum numbers identical to the squarks d and u. 
For these leptoquarks the couplings A correspond to the Yukawa couplings, \[-f,, which describe 
the I^p supersymmetric LiQjDk interaction [18]. Here i, j and k are the family indices and 
Li, Qj and D^ are the super-fields containing the left-handed leptons, the left-handed up-type 
quarks and the right-handed down-type quarks, respectively, together with their supersymmetric 
partners. 



2.3 Large extra dimensions 

In some string inspired models the small nature of the gravitational force is explained by the 
existence of compactified extra dimensions |fT9l . In these models the gravitational scale Ms in 
4 + ra dimensions is related to the size R of the compactified extra dimensions via the Planck 
scale Mp ~ R"- Ml^"". SM particles reside on a four-dimensional brane, while the spin 2 
graviton propagates into the extra spatial dimensions creating a tower of Kaluza- Klein states. 
Assuming that the ultraviolet cut-off scale of the tower is of similar size to the gravitational 
scale, an effective contact-type interaction [|20ll term can be defined with a coupling coefficient: 

The coupling A depends on details of the theory and is conventionally set to ±1. 



2.4 Quark Radius 

A clear manifestation of substructure would be the observation of finite size effects like the 
measurement of electroweak charge distributions of fermions. Finite size effects are typically 
described by a standard form factor in the eg scattering cross section: 

f{Q') = l-^Q^ (5) 

o 

which relates the decrease of the scattering cross sections at high Q^ to the mean squared radius 
(i?^) of the electroweak charge distribution. This form factor modifies the Q^ dependence of 
the ep scattering cross section similarly to the CI models described above. 



3 Data and Analysis Method 

The analysed data sample is recorded in e^p and e^p collisions corresponding to an integrated 
luminosity of 281 pb^^ and 165 pb^^, respectively. The measurement of the differential neutral 
current cross-section, da/dQ"^ , which is used to probe possible CI signatures follows the pre- 
vious measurements based on data recorded in the years 1994-2000 [!2Tl - |23l and includes new 
data recorded from 2003-2007. A list of the analysed data sets is given in tablelB 



Reaction 


Cint [pb ^ 


v^ [GeV] 


Polarisation (Pe [%]) 


e^p — )■ e+X 


36 


301 


Unpolarised 


e~p — )• e~X 


16 


319 


Unpolarised 


e^p — )• e+X 


65 


319 


Unpolarised 


e~p —7- e^X 


46 


319 


Right (Pe = +37) 


e~p — )■ e~X 


103 


319 


Left (Pe = -26) 


e^p — )■ e+X 


98 


319 


Right (P, = +33) 


e+p -)■ e+X 


82 


319 


Left (Pe = -38) 



Table 1: Data samples recorded in the years 1994-2007 with corresponding integrated lumi- 
nosities, centre-of-mass energies and average longitudinal polarisations. 

The data collected from the year 2003 onwards were taken with a longitudinally polarised 
lepton beam, with typical polarisation values of ±35%. The average luminosity weighted po- 
larisations of the e^p and e^p data sets are small. Due to the different chiral structure of hypo- 
thetical new particles with respect to the irreducible SM background, the sensitivity to possible 
new physics phenomena is increased by up to 15 percent by analysing the data sets with left and 
right longitudinal lepton polarisation separately. 

Contact interactions are investigated by searching for deviations in the NC differential cross 
section da/dQ"^ from the SM expectation at high negative four-momentum transfer squared 
Q^ > 200 GeV^. The SM cross section of neutral current scattering factorises into the elec- 
troweak matrix element of the hard eg interaction process and the parton distribution function 
(PDF) of the proton. The Q'^ dependence of the PDF is calculated using perturbative QCD [|24ll . 



For the CI analysis the parton densities at high values of Q^, corresponding to high values of 
X, are of special importance. In this analysis the CTEQ6m [25 1 PDF is used to calculate both the 
SM and signal expectations. The CTEQ6m set was obtained by fitting several experimental data 
sets. At high x this PDF is mostly constrained by fixed target experiments and also by VT-boson 
production and jet data from the Tevatron experiments, which are not sensitive to possible eq 
contact interaction processes. CTEQ6m also includes early e^p scattering data at high Q'^ from 
the HI {C = 52 pb^^) and ZEUS (C = 30 pb^^) experiments. However, since the e+p (e~p) 
data sets analysed here are 6(10) times larger, the residual correlations between the HERA data 
and the CTEQ6m PDF are small and are neglected in the following. Furthermore, the CTEQ6m 
parton densities can be regarded as unbiased with respect to possible contact interaction effects. 
CTEQ6m is chosen as it describes many experimental data and in particular, the HERA data in 
the region Q^ < 200 GeV^, which are not used in this analysis. The results of this analysis are 
verified using an alternative PDF not based on HERA high Q'^ data, as described in section |4l 

The single differential cross sections da/dQ'^ are measured for e^p and e^p scattering up 
to Q"^ = 30000 GeV^ and compared to the SM expectation. The ratio of the data to the SM 
expectation is shown in figure [TJ Good agreement between data and the SM is observed, in 
particular in the high Q^ region, which is the focus of this analysis. 

In the next step, a quantitative test of the SM and the CI models is performed by investigating 
the measured cross sections da/dQ"^ following the analysis method described in [|T1 by applying 
a minimisation of the x^ function [|26ll : 



.. . _ V (^r - cr^v) (1 - E. A..(..)))^ 



^Ltat (^r^^l^iv) (1 - Efe ^ikiSk)) + (^i.uncor ^ ■ 



ei . (6) 



k 



Here a'^^^ and af'^i]) are the experimental and theoretical cross sections, respectively, for the 
measurement point i, and i^j stat and 5j uncor correspond to the relative statistical and uncorrelated 
systematic errors, respectively. The theoretical cross section includes both the SM and the con- 
tact interaction term, and it depends on the coupling coefficient i], which is varied in the fit. The 
functions Aik^Sk) describe the correlated systematic errors for point i associated to a source k 
and depend on the fit parameters Sk- In particular, the normalisations of the individual data sets 
described in table[T]are free parameters, only constrained by the individual luminosity measure- 
ments. Since the precise cross section measurements at low Q^ determine the normalisations, 
new physics signals are mainly tested by exploiting the shape of the Q^ distribution. 

Statistical and systematic uncertainties are taken into account in the fit procedure. The 
following sources of experimental uncertainties are accounted for [|2TI - l23ll27l : the electromag- 
netic energy scale uncertainty of 1 — 3%, the polar angle uncertainty of the scattered lepton of 
2 — 3 mrad, the uncertainty on the electron identification of 0.5 — 2%, the hadronic energy scale 
uncertainty of 2 — 7% the uncertainty from the luminosity measurement of 1.6 — 3.8% and the 
uncertainty on the electron beam polarisation of 1 — 2.3%. The effect of the above systematic 
uncertainties on the SM expectation is determined by varying the experimental quantities by 
±1 standard deviation in the MC samples and propagating these variations through the whole 
analysis. The uncorrelated systematic uncertainties of the measurements vary as function ofQ^ 
between 1 — 11% (1.6 — 13%) for e~^p {e~p) scattering. The dominant sources of the correlated 
systematic errors are the PDF uncertainty (about 8%), the uncertainties from the luminosity 
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measurement and from the experimental uncertainties on the energy scale and the polar angle 
of the scattered lepton. All other experimental systematic uncertainties are found to have a 
negligible impact on the analysis. 



4 Results 

The data used in this analysis are found to be consistent with the expectation from the SM alone 
(r/ = in equation[6l) based on the CTEQ6m PDF, yielding a x^/dof = 16.4/17 (7.0/17) for 
the e^p (e~p) data. The normalisation constants of the individual data sets agree well with the 
SM expectation within the PDF uncertainties. 

For each CI model the effective scale parameters and couplings describing the new physics 
scale are determined by a fit to the differential NC cross section. All scale parameters are found 
to be consistent with the SM and limits are calculated at 95% confidence levels (CL) using the 
frequentist method as described in the previous publication ^. 

Lower limits on the compositeness scale A in the context of the general contact interaction 
model are presented in table [2] and figure [21 The results are presented for eight scenarios, which 
differ in their chiral structure as determined by the CI coupling coefficients r^^^. Depending 
on the model and the sign of the coefficients, limits on A are obtained in the range 3.6 TeV 
to 7.2 TeV. In figure |3l differential cross section measurements for e^p and e^p scattering 
normalised to the SM expectation are compared to the predictions corresponding to the 95% 
CL exclusion limits of the VV model, A^^ > 5.6 TeV and Ayy > 7.2 TeV. 

For leptoquark-type contact interactions, the notation, quantum numbers and lower limits on 
Mlq/A are presented in table[3l The limits are in the range Mlq/A > 0.41 — 1.86 TeV. Lepto- 
quarks coupling to u quarks are probed with higher sensitivity, corresponding to more stringent 
limits than those coupling to d quarks due to the different quark densities in the proton. In fig- 
ureHl the normalised differential e'^p cross section measurements are compared to the predicted 
cross sections corresponding to the 95% CL exclusion limits of the S^ and V-^ leptoquarks. At 
high Q^ the existence of a S'f (Vl") leptoquark would lead to an increase (decrease) of the 
e^p cross sections, which is not observed. For a Yukawa coupling of electromagnetic strength, 
A = 0.3, scalar and vector leptoquark masses up to 0.33 TeV and 0.56 TeV are excluded, re- 
spectively, comparable or exceeding limits obtained by Tevatron and LHC. The leptoquarks Sq 
and S^,2 may also be interpreted as squarks in the framework of ^p supersymmetry and the 
corresponding limits in terms of the ratio M^/X' are given in table IH For a Yukawa coupling of 
electromagnetic strength, the corresponding lower limit on the u mass of 0.33 TeV is similar to 
that obtained recently by the HI collaboration in a direct search Q. 

Lower limits in a model with large extra dimensions on the gravitational scale Ms in A+n di- 
mensions assuming a positive (A = +1) or negative (A = —1) coupling are given in tabled 
Mass scales Ms < 0.9 GeV are excluded at 95% CL. The corresponding cross section predic- 
tions normalised to the SM expectation are compared to the e^p data in figure [5l 

Finally, an upper limit at 95% CL on the quark radius R^ < 0.65 ■ 10"^*^ m is derived 
assuming point-like leptons. The corresponding cross section predictions normalised to the SM 
expectation are compared to the e^p data in figure [6l 
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The above results are also verified using a dedicated HI PDF set based on data collected in 
the years 1994-2007. This PDF set was obtained from a next-to-leading order QCD fit to the 
HI data [26] with Q^ < 200 GeV^, excluding the high Q'^ data used in this analysis. Both the 
SM expectation and limits derived using the dedicated HI PDF agree well with those obtained 
using the CTEQ6m PDF within the uncertainties. 



5 Summary 

Neutral current deep inelastic e~p and e+p scattering cross section measurements are analysed 
to search for new phenomena mediated via contact interactions. The data are well described 
by the Standard Model expectations. Limits on the parameters of various contact interaction 
models are presented at 95% CL. 

Lower limits on the compositeness scale A are derived within a general contact interaction 
analysis. The limits range between 3.6 TeV and 7.2 TeV depending on the chiral structure, 
corresponding to an increase by a factor of about two compared to previous HERA searches. 
The study of leptoquark exchange yields lower limits on the ratio Mlq/A between 0.41 TeV 
and L86 TeV, considerably improving constraints from the previous analysis. Squarks in the 
framework of -R-parity violating supersymmetry with masses satisfying Mu/\[-^ < 1.10 TeV 
and M^/\\^^ < 0.66 TeV are excluded. Possible effects of low-scale quantum gravity with 
gravitons propagating into extra spatial dimensions are also investigated, where lower limits on 
the gravitational scale in 4 + n dimensions Ms > 0.9 TeV are found. Finally, a form factor 
approach yields an upper limit on the size of light u and d quarks of Rg < 0.65 • 10~^^ m, 
assuming point-like leptons. 

Using the full HERA data set, limits derived in this analysis are more stringent than previous 
results by HI and ZEUS. The results can also be compared to those obtained by the LEP, 
Tevatron and, most recently, LHC collaborations. For most models with Yukawa couplings of 
electromagnetic strength, or stronger, the analysis presented here provides the most stringent 
limits on first generation leptoquarks. 
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HI Search for General Compositeness 
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Table 2: Lower limits at 95% CL on the compositeness scale A. The A+ limits correspond 
to the upper signs and the A^ limits correspond to the lower signs of the chiral coefficients 

l-^LL^ e-LR^ e-RL^ ^RRi- 
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HI Search for Heavy Leptoquarks 
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Table 3: Lower limits at 95% CL on Mlq/A for scalar (S) and vector (V) leptoquarks, where 
L and R denote the lepton chirality and the subscript (0, 1/2, 1) is the weak isospin. For each 
leptoquark type, the relevant coefficients e^^ and fermion number F = L + 3B are indicated. 
Leptoquarks with identical quantum numbers except for weak hypercharge are distinguished 
using a tilde, for example Vq^ and Vq^. Quantum numbers and helicities refer to e^q and e^q 
states. 



HI Search for^ Squarks 

Channel Coupling e^^ Af,-/A' [TeV] 

e+c^^iiW All, e^^ = +i 1.10 

e-u -> (i(^) a;,! 4r = -1 0.66 



Table 4: Lower limits at 95% CL on Mq/\' for the Rp violating couplings X'^^^ |T8l, where i, j, k 
are family indices. The coefficients e^j^ are also shown. The A'^^^^ i^'iji) coupling corresponds to 
the Slj' {S^,^) leptoquark coupling shown in tabled 
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HI Search for Large Extra Dimensions 

VG = X/Ml 
coupling A Ms [TeV] 

+1 0.90 

-1 0.92 



Table 5: Lower limits at 95% CL on a model with large extra dimensions on the gravitational 
scale Ms in 4 + n dimensions, assuming positive (A = +1) or negative (A = —1) couplings. 
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Figure 1 : The ratio of the measured cross section to the Standard Model prediction determined 
using the CTEQ6m PDF set for e+p — )■ e+X and e~p — > e^X scattering. The top figurecor- 
responds to the full HI data with an average longitudinal polarisation of P ^ 0. The middle 
and bottom figuresrepresent polarised HI data taken from the year 2003 onwards for different 
lepton charge and polarisation data sets. The error bars represent the statistical and uncorrelated 
systematic errors added in quadrature. The bands indicate the PDF uncertainties of the Standard 
Model cross section predictions. 
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H1 Search for General Compositeness 




A- [TeV] 



Figure 2: Lower limits at 95% CL on the compositeness scale A for various chiral models, 
obtained from the full HI data. Limits are given for both signs A"*" and A~ of the chiral coeffi- 
cients. 
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Search for General Compositeness 
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Figure 3: The measured neutral current cross section dcr/dQ^ normalised to the Standard Model 
expectation. HI e'^p scattering data are compared with curves corresponding to 95% CL exclu- 
sion limits obtained from the full HI data for the VV compositeness scale model, for both signs 
A+ and A^ of the chiral coefficients. The error bars represent the statistical and uncorrected 
systematic errors added in quadrature. 
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Figure 4: The measured neutral current cross section dcr/dQ^ normalised to the Standard Model 
expectation. HI e^p scattering data are compared with curves corresponding to 95% CL exclu- 
sion limits obtained from the full HI data on the ratio Mlq/A for the S\ and V/" leptoquarks. 
The error bars represent the statistical and uncorrected systematic errors added in quadrature. 
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Search for Large Extra Dimensions 
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Figure 5: The measured neutral current cross section dcr/dQ^ normalised to the Standard Model 
expectation. HI e^p scattering data are compared with curves corresponding to 95% CL ex- 
clusion limits obtained from the full HI data on the gravitational scale, Ms for both positive 
(A = +1) and negative (A = —1) couplings. The error bars represent the statistical and uncor- 
related systematic errors added in quadrature. 
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Figure 6: The measured neutral current cross section da/dQ'^ normalised to the Standard Model 
expectation. HI e^p scattering data are compared with curves corresponding to 95% CL exclu- 
sion limits obtained from the full HI data on the quark radius, Rg assuming point-like leptons. 
The error bars represent the statistical and uncorrected systematic errors added in quadrature. 
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